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I discuss a study performed by the SCRI lattice gauge theory group which compared light hadron spectroscopy 
using the Wilson and Clover fermionic actions. We showed that a Clover coefficient chosen using tadpole-improved 
tree-level perturbation theory effectively eliminates the 0{a) discretization errors present in the Wilson action. 
We found that discretization errors in light spectroscopy for both the Wilson and Clover actions are characterized 
by an energy scale of about 200-300 MeV, indicating that these errors can be reduced to the 5% level by using 
the Clover action at an inverse lattice spacing of about 1.3 GeV. 



1. INTRODUCTION 

In this talk, I will discuss the results of an ongo- 
ing project at SCRI to understand the systematic 
errors in lattice light spectroscopy calculations. 

The short term goal of our program has been 
to determine the value of the lattice spacing at 
which the Clover lattice fermion action reduces 
discretization errors to the 5% level. Once this 
lattice spacing is known, we can perform dynam- 
ical simulations on a large physical volume and 
attempt to match our results to experiment. 

To determine this lattice spacing, we have con- 
centrated on two sub-projects. First, we have 
performed a detailed comparison of light spec- 
troscopy using the Wilson and Clover discretiza- 
tions on a pre-existing gluonic ensemble, gener- 
ated by the HEMCGC collaboration, which has 
an inverse lattice spacing of about 2 GeV and 
includes the effects of two flavors of (staggered) 
light quarks. The purpose of this study was to 
see if, at typical lattice spacings, the Clover dis- 
cretization significantly improved the discrepan- 
cies with experiment found in Wilson simulations, 
particularly in ratios of observables which probe 
different energy scales (such as the ratio of bot- 
tomonium splittings to the proton mass). Sec- 
ond, we have used an improved gluonic action on 
coarse lattices to determine the size of the scaling 
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violations directly where they are large. 

This talk is organized as follows. In section 
2, I discuss technical details of our calculations. 
In section 3, I discuss the techniques we used in 
the scaling analysis of our quenched spectroscopy 
data, while in section 4, I present the results of 
these calculations. In section 5, I present the 
results of our comparison of Wilson and Clover 
spectroscopy, including the effects of two flavors 
of dynamical fermions, at an intermediate lattice 
spacing. In section 6, I summarize our conclu- 
sions. 

2. THE SIMULATIONS AND FITTING 
PROCEDURES 

We used two different types of gluonic ensemble 
in this study. To determine the size of discretiza- 
tion effects, we generated quenched ensembles us- 
ing a one loop 0(a 2 ) improved gluonic action 
at six values of /3. The inverse lattice spacing of 
these ensembles ranged from 1.3 to 0.5 GeV, cor- 
responding to lattice spacings between .15 and 0.4 
fermi (see table [l]). The ensemble at each j3 con- 
sisted of 100 configurations of size 16 3 x 32. For 
the dynamical runs, we used 16 3 x 32 gauge con- 
figurations generated by the HEMCGC collabo- 
ration ||, with (3 = 5.6 and 2 flavors of staggered 
fermions at a mass of am = 0.01. The inverse lat- 
tice spacing for this ensemble has previously been 
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determined to be roughly 2 GeV. 



Table 1 

Values of 0(a 2 ) improved f3 and corresponding 
inverse lattice spacings for quenched study. 






7.9 


7.75 


7.6 


7.4 


7.1 


6.8 


a" 1 (GeV) 


1.3 


1.1 


.96 


.78 


.59 


.48 



On all ensembles, both quenched and dynam- 
ical, we calculated light hadron correlation func- 
tions using both Wilson and tadpole improved 
Clover quark propagators. In all cases we used 
methods involving multi-state fits to multiple cor- 
relation functions to extract hadron masses, the 
technical details of which are discussed in H . 

The main focus of this work, especially of 
the quenched portion, was to determine the size 
of scaling violation effects rather than to at- 
tempt to reproduce or predict experimental re- 
sults. Since discretization errors should get worse 
with increasing quark mass (because the typical 
energy scales should grow), we did not devote 
much effort to chiral extrapolations. Most of the 
quenched results I will present here used a low- 
est order (quadratic in Mps) chiral fit ansatz, 
while the dynamical results were analyzed using 
both quadratic and cubic ansatze. Gottlieb has 
already presented a detailed discussion of this ex- 
tremely important (and ill-controlled) source of 
systematic errors in his talk ||, so I will not dis- 
cuss chiral extrapolations here, other than to say 
that our results agree qualitatively with his. 

2.1. Low Quark Masses on Coarse Lattices 

I would like to make two technical comments 
concerning coarse lattices before proceeding to 
the bulk of the talk. First, for a particular di- 
mensionful value of Mps, going to a coarse lat- 
tice does not seem to improve the speed of con- 
vergence of the fermion inversion algorithm. In 
other words, the condition number of the fermion 
kernel seems to be controlled by m q j Aqcd, the 
quark mass over the QCD scale, rather than am q , 
the dimensionless quark mass. This is a very rea- 



sonable behavior, but it means that the gains to 
be found in going to a coarse lattice might be 
smaller than one might think. 

The other comment involves exceptional config- 
urations. We had originally planned to perform 
our quenched coarse lattice runs using 8 3 x 16 con- 
figurations. At low values of /3, however, we found 
many exceptional configurations. This forced us 
to increase the lattice size to 16 3 x 32, which 
suppressed the effect. This method of increasing 
the volume to reduce exceptional configurations, 
can, however, result in an unacceptable increase 
in computational cost (which scales linearly with 
the lattice volume in the regime where finite vol- 
ume effects are unimportant) unless one is care- 
ful. This is because each inversion of a local- 
ized hadronic source will only sample a roughly 
(lfm) 3 volume of the lattice, so there is no gain in 
statistics when one increases the lattice volume. 
The solution is to smear the hadron around mul- 
tiple origins (separated by some multiple of the 
hadron's radius) on the initial timcslice. Naively 
this would increase the computational cost by 
the number of origins used, but one can invert 
the quark propagators from all origins simultane- 
ously by superposing them with a random Z(3) 
phase ||. When using this method the compu- 
tational cost of inverting from multiple origins is 
the same as the cost to invert from a single origin. 

In practice, this means that one can double 
the spatial extent of the lattice (decreasing both 
finite volume and exceptional configuration ef- 
fects) without increasing the computational cost 
of achieving equal statisticsl This is because, al- 
though each inversion will be eight times as ex- 
pensive, it will also yield eight times as many sta- 
tistically independent measurements. This means 
that (ignoring memory limitations, thermaliza- 
tion times, and the reduction of HMC acceptance 
rates on larger volumes) one can always effectively 
eliminate finite volume effects from one's simula- 
tions at no additional CPU cost by increasing the 
lattice volume and using the superposition trick. 
Although the superposition trick is vital on coarse 
lattices (where the lattice volume must be large), 
it should also be used to eliminate finite volume 
effects on fine lattices. 
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3. SCALING FITS 

The main idea in performing extrapolations to 
the continuum limit, as well as in attempting 
to reduce discretization errors through improve- 
ment programs, is to Taylor expand dimensionless 
physical ratios around a — 0. For example, 
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=M(0) 



M'(0), „ M"(0). „j 



:M(0) fl + C (1) a + C (2) a 2 + 
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where M(0), Af'(0), . . . are the continuum limits 
of the ratio and its derivatives (note that I will 
work in units of the string tension a when dis- 
cussing scaling, even though it is not an exper- 
imental observable). Continuum extrapolations 
attempt to truncate the expansion at some order 
(usually first) and fit to data obtained at different 
a's to determine M(0). Improvement programs 
modify the action so as to reduce or eliminate the 
low order coefficients such as C' 1 ' and C^ 2 \ For 
the (non-perturbatively tuned) Clover action, the 
leading discretization error is C< 2 V. Our strat- 
egy in this project has been to determine by 
running on coarse lattices with large a. Once we 
have determined C' 2 ' , we can choose a value of a 
for our production runs at which the discretiza- 
tion errors are comparable to the other errors in 
our simulation. In particular, it is a waste of 
computational effort to run at values of a where 
the attainable statistical errors are an order of 
magnitude larger than the discretization errors; 
higher statistics can be obtained for the same ef- 
fort on coarser lattices, while the discretization 
errors will still be negligible. 

3.1. fi 

Discretization errors in lattice calculations are 
caused by irrelevant, higher dimension opera- 
tors appearing in the cut-off action. In general, 
physical observables will have some characteristic 
infra-red energy scale, fi, which is the typical en- 
ergy scale of the quantum fluctuations important 
to the observable. For light hadronic observables, 
we expect this scale to be set by Aqcd or the 



constituent quark mass, i.e. /i should be several 
hundred MeV. 

If we rewrite the dimensionful coefficients of 
equation [l], C, in terms of dimensionless coeffi- 
cients D, i.e. C(") = D^fi n : then a reasonable 
discretization should have D's which are 0(1). If 
we assume that this is the case (i.e. set D = 1), 
then each coefficient in equation [P gives us an in- 
dependent estimate of /i: 



/'r. 



n)\l/n 



(2) 



If the different estimates /i n are in rough agree- 
ment, then it is reasonable to estimate the higher 
order 0(a N ) corrections to be about (n n a) N . 
Conversely, if the estimates \x n wildly disagree, 
then we know that the coefficients of the a ex- 
pansion are not very well behaved, and we should 
be very cautious in inferring the size of higher 
order (unmeasured) discretization errors. In our 
quenched study, we found that /i ra 's for the p, the 
nuclcon, and the A were all in the region 200-300 
MeV, for both the Wilson and Clover discretiza- 
tions. This means that one can expect discretiza- 
tion errors of about 1-2% when using the Clover 
action at a lattice spacing of .1 fm. If one wants 
discretization errors to be about 5%, as is reason- 
able for a first effort at reproducing the physical 
spectrum using dynamical fermions, then a lat- 
tice spacing of about .2 fm appears to be the best 
choice. 

Note that [i is closely related to the perturba- 
tive q* of Lepage and Mackenzie |Q. In particular, 
one might consider equation^, for a particular n, 
to be a non-perturbative definition of q*. It is 
not necessary, however, that the physical scale fx 
be close to our estimates /x„; a particularly bad 
discretization of the action, for example, would 
yield discretization errors much worse than the 
expected (/ia)", leading to spuriously large val- 
ues of For example, one might expect the 
/x n 's in calculations using Staggered fermions to 
be about twice those of Wilson or Clover calcula- 
tions because of the two hop nature of the Stag- 
gered derivative operator. It is likely, however, 
that fi is a (rough) lower bound on /x„, since it is 
hard to imagine an action with spuriously small 
values of the D's. 
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3.2. Scaling Fit Ansatze 

We used linear and quadratic ansatze for our 
scaling fits of Wilson and Clover data. We called 
these ansatze 1, 2, and 12, corresponding to pure 
linear (1), pure quadratic (2), or mixed linear plus 
quadratic (12). All of our scaling fits are done in 
units of the string tension; for physical predic- 
tions we must take ratios of extrapolated values 
(to remove the string tension, which is not exper- 
imentally measurable). 

The first question we attempted to answer was 
whether the Clover action with a tadpole im- 
proved tree level coefficient is effective in re- 
moving the linear discretization errors of Wilson 
fermions. If the Clover coefficient, Csw is non- 
perturbatively tuned, then the lowest order dis- 
cretization errors will be quadratic. In our calcu- 
lations, however, we used tadpole-improved tree 
level perturbation theory to set Csw, so we ex- 
pect 0(a(a)a) errors in addition to 0(a 2 ). Since 
a(a) runs fairly slowly, we treated it as a constant 
over the range of lattice spacings in our calcula- 
tion, i.e. we included it in our fit ansatze as a 
linear term with a coefficient which we expect to 
be small (0(a)). 



Table 2 

Q values for linear, quadratic and mixed scaling 
fits to Wilson and Clover vector meson masses. 
V indicates the benchmark vector particle at 
Mp S /M v = 0.7. 



State 


Action 


linear 


quadratic 


mixed 


P 


Wilson 


.41 


« lO" 6 


.34 


Clover 


.20 


.98 


.94 


V 


Wilson 


.28 


« lO- 1 " 


.45 


Clover 


.04 


.94 


.94 



Table g contains the values of the confidence 
level, Q, for scaling fits of the vector mass at all 
six (3 values, using the 1, 2, and 12 ansatze. Both 
Wilson and Clover results are shown for vector 
masses corresponding to the p (i.e. Mps/My = 
Mir/Mp) and to a benchmark value agreed to at 



the recent Seattle workshop g, for which chi- 
ral extrapolations are not required (Mps/My = 
0.7). For the Wilson fits it is clear, as expected, 
that a linear term is necessary; the pure quadratic 
ansatz is clearly ruled out by the small Q values. 
For the Clover fits, however, one cannot draw a 
definitive conclusion from these Q values alone; 
both the linear and quadratic fits to the p have 
acceptable Q, although the linear Q's (especially 
.04 for the V) are low. The message to be drawn 
from Table || is that it is extremely difficult to 
distinguish between analytic functional forms to 
be used in extrapolations by looking at confidence 
levels for fits to a single curve; roughly 2 Gflops-yr 
of CPU time were required to generate the Clover 
data yet these fits still are unable to distinguish 
linear from quadratic behavior. This observation 
is important not only when doing lattice spacing 
extrapolations, but also when doing chiral extrap- 
olations. I will show later that simultaneous fits 
to both Wilson and Clover data do rule out the 
pure linear Clover ansatz; unfortunately we do 
not know how to apply this technique to chiral 
extrapolations. 

Although the linear and quadratic ansatze can- 
not be distinguished by Q, they result in signif- 
icantly different continuum extrapolations. Fig- 
ure [j] shows the results of the three different fits 
to the Clover M p data, along with the contin- 
uum extrapolations. It is obvious that the linear 
and quadratic extrapolations differ by many stan- 
dard devations. This is a serious problem, since 
it shows that there is a systematic error associ- 
ated with choice of fit ansatz which can be many 
(statistical) a. 

Our analytic knowledge tells us that the (per- 
turbatively tuned) Clover action should reduce 
but not eliminate the 0(a) errors. The mixed 
ansatz permits both the large quadratric and sup- 
pressed linear terms to appear in the extrapola- 
tion. As can be seen from figure |l|, the linear coef- 
ficient is consistent with zero but the extrapolated 
value has a much larger statistical uncertainty. 
This is because (as in all other cases in table ||) 
the mixed fit is underdetermined; an ansatz with 
fewer parameters (linear for Wilson or quadratic 
for Clover) fits the same data with good Q. We 
expect the true uncertainty of the Clover extrap- 



5 



3.0 



1 6 



1.4 



1.2 



■ i \^ 




-i - :-: 


x Clover — 











0.0 



0.5 



1.0 



Table 3 

Q values for joint scaling fits to Wilson and Clover 
vector meson masses. V indicates the benchmark 
vector particle at MpsjMy = 0.7. 



State 


Ansatz 


Wilson 


Clover 


linear 


quadratic 


mixed 


9 


linear 


« 10~ 2ti 


.016 


.69 


mixed 


.002 


.85 


.78 


V 


linear 


« io- ,/a 


« 10-^ 


.16 


mixed 


.001 


.98 


.96 



Figure 1. Results of different scaling fits to the 
p mass. 



olation to fall somewhere between the quadratic 
and mixed error bars. In presenting extrapolated 
values I will quote both estimates of the uncer- 
tainty. 

3.3. Joint Fits 

Until now, I have presented the Clover and Wil- 
son results separately. The zeroth order coeffi- 
cient in both cases, however, is just the contin- 
uum limit of the quantity being fit. This value 
should not depend on the discretization scheme 
used, so we can obtain more constrained fits by 
simultaneously fitting Clover and Wilson data: 

M w (a^) = M(0) [1 + Oufa) + {»fa) 2 ] 

Mc(aV^) = M(0) [1 + (p%a) 2 ] , (3) 

where M is the quantity being fit (i.e. aMy /ay/a) 
and the subscript indicates data from Wilson (W) 
or Clover (C) discretizations. In this example, I 
have used the 2 ansatz for the Clover data, and 
the 12 ansatz for Wilson. Although the Wilson 
and Clover data at a particular value of a are cor- 
related, we have not included these correlations in 
our joint fits. This could lead us, in principle, to 
overestimate the confidence levels of our fits, but 
we expect this effect to be small. 



Table || presents the Q's for joint fits using all 
three ansatze for Clover and the linear and mixed 
ansatze for Wilson data. The quadratic Wilson 
ansatz is not considered because it has already 
been ruled out in table g. It is now obvious that 
the pure linear ansatz for Clover is also ruled out; 
its extrapolated value cannot be made to agree 
with that of any of the Wilson ansatze. Although 
the 12_1 (mixed Clover, linear Wilson) fits have 
reasonable Q, they are unstable when more pa- 
rameters are added; the parameters differ from 
those in the 12.12 fits by 1 — 2 standard devia- 
tions for the p and 3a for the V (using the under- 
determined 12_12 error bars). As we saw when 
comparing the Clover fits in table |^, this is an 
indication of an incorrect ansatz. This behavior 
should be compared with that of the 2_12 fits, all 
of whose parameters agree with those in the 12.12 
fits. This is entirely consistent with our analytic 
expectations; we are at a coarse enough lattice 
spacing that we can resolve the quadratic errors 
in Wilson fermions in addition to the leading lin- 
ear errors, which the Clover term suppressed suf- 
ficiently to allow a pure quadratic ansatz. When 
an additional linear term is included (i.e. 12.12 
rather than 2.12 ansatz), the parameters are es- 
sentially unchanged, albeit with larger error bars, 
while the linear term is consistent with zero. 

In presenting the scaling results that follow, I 
will use central values from the 2_12 fits and quote 
error bars from both the 2.12 and 12.12 fits as 
"best case" and "worst case" uncertainties. 
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4. SCALING RESULTS 



cl 2, wl 12 




0.0 0.5 1.0 

(7a 

Figure 2. Joint scaling fit to the benchmark 
vector mass, using a pure quadratic ansatz for 
Clover(+) and mixed ansatz for Wilson(o) data. 



Figure || shows the results of a joint 2_12 fit to 
the benchmark vector mass, with the string ten- 
sion used to set the scale. The most striking fea- 
ture is that, at reasonable (3, the Clover term has 
succeeded in eliminating most of the discretiza- 
tion effects, even though its coefficient was set 
using tree-level (tadpole-improved) perturbation 
theory. Note that the two highest (3 values cor- 
respond to a inverse lattice spacing of about f 
GeV; a 6 GeV inverse lattice spacing would be 
required before the Wilson discretization errors 
are reduced to the same level. The slope param- 
eters from the fit, p% , pY and pY as defined in 
equation [| are 240, 300, and 180 MeV, respec- 
tively, where we have taken y/a to be 440 MeV. 
The fact that all three estimates are about the 
same size is an indication that the discretization 
scale p is roughly 200-300MeV. 

We have performed similar fits to the p mass, 
as well as to the masses of the nucleon and the 
A and the nucleon-p mass ratio. The results of 
these fits are presented in table ||, and seem to 
present a consistent picture. The V corresponds 



to a quark mass near the strange, so one expects 
the //'s to be harder than for the p, which con- 
tains massless quarks. The p's in table || of the 
V are all slightly harder than those of the p. Al- 
though this difference is about the same size as 
the uncertainty in the p estimates, the different 
vector results are correlated so the effect is prob- 
ably significant. The baryonic /i's, in the range 
150-200MeV, seem to be slightly smaller than 
those of the vectors, which would lead to smaller 
discretization errors at a particular lattice spac- 
ing. This is born out by our dynamical fermion 
study (discussed below) , where we were unable to 
distinguish between Clover and Wilson baryonic 
masses, but were able to see a difference in the 
p mass. The p estimate coming from the linear 
term in the Wilson nucleon-p mass ratio at first 
appears small, but this is due to the fact that 
the linear coefficient for each particle is negative, 
leading to a large cancellation in the ratio. 

Table 4 

Estimates of p from coefficients of joint scaling 
fits to Wilson and Clover results for various ob- 
servables. Masses are measured in units of the 
square root of the string tension, which was taken 
to be 440 MeV when estimating p. V indicates 
the benchmark vector particle at Mps /My = 0.7. 



Observable 


Clover 


Wilson 




n w 

n 


V$ 


M p 


220 


250 


130 


My 


240 


300 


180 


M N 


160 


190 


160 


Ma 


190 


190 


80 


M N /M P 


180 


50 


230 



Overall, we conclude that the p for light 
hadron spectroscopy, using Wilson- type actions, 
is roughly 200-300MeV. If we take a worst-case 
estimate of 300MeV and assume pure quadratic 
errors for the Clover action, this implies that lat- 
tice spacings of 2.1 and 1.3 GeV are required to 
obtain discretization errors of 2% and 5%, respec- 
tively. To obtain the same errors with the Wil- 
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son action, one would need corresponding lattice 
spacings of 15 and 6 GeV, respectively! 




0.0 0.2 0.4 



Figure 3. World scaling data for the p mass. 
IG indicates an improved gluonic action (our re- 
sults), while NP uses the chiral ward identity to 
tune the Clover parameter. 



To illustrate the effectiveness of using the 
Clover action at intermediate lattice spacings, fig- 
ure H shows a sample of world data for the p mass. 
The discretization errors in our Clover point at 
(3 = 7.9 (corresponding to an inverse lattice spac- 
ing of about 1.3 GeV) are about the same size as 
the statistical errors, and are comparable to the 
statistical errors of other calculations. 

Although it was not the main goal of this 
project, we have obtained the quenched contin- 
uum values of the p, nucleon, A, and 0.7 bench- 
mark vector masses in units of y/a, as well as the 
nucleon-p mass ratio from our joint fits. These 
results are presented in table |^, with best and 
worst case error bars shown. The benchmark 
V was chosen to be in the vicinity of computa- 
tionally accessible quark masses; it is unlikely to 
be sensitive to chiral extrapolation uncertainties. 
The physical hadrons, however, required a large 



Table 5 

Continuum values for quenched light hadron 
masses and the nucleon-p mass ratio. Masses are 
given in units of the string tension. Statistical er- 
rors from both 2_12 and 12_12 extrapolations are 
shown. V indicates the benchmark vector particle 
at Mps/My = 0.7. Systematic chiral extrapola- 
tion errors for physical masses are not estimated 
(for My they should be negligible). 



Observable 


Value (y/cr) 


M v 


2.122(11)(56) 


M p 


1.758(11)(46) 


M N 


2.30(3)(12) 


M A 


2.93(5)(15) 


M N /M p 


1.297(17)(81) 



extrapolation (using an ansatz linear in Mp s )\ 
these results are susceptible to (undetermined) 
systematic errors arising from other terms in the 
chiral expansion. The only value which can be 
compared directly to experiment is the nucleon- 
p mass ratio. We obtain 1.297(17)(81), which is 
consistent with the experimental value of 1.22 if 
one uses the worst case uncertainty. 

5. DYNAMICAL RESULTS 

Our dynamical runs were intended to investi- 
gate whether Clover results on a fine lattice with 
dynamical fermions (partially) included would re- 
solve discrepancies with the experimental spec- 
trum found in Wilson studies, and to check how 
large the differences between Wilson and Clover 
spectroscopy are on a fine lattice. The most im- 
portant of these discrepancies is the disagreement 
between the lattice spacings obtained from light 
spectroscopy and those obtained from T split- 
tings |J. Previous Wilson calculations on this 
ensemble jl0| found inverse lattice spacings of 
about 2 GeV. Spin-independent bottomonium 
splittings calculated on the same ensemble using 
the NRQCD formalism, however, obtained values 
near 2.4 GeV - a 20% discrepancy. It is unlikely 
that this would change if three flavors of dynam- 
ical fermion were used, since quenched results at 



8 



a similar lattice spacing show roughly the same 
discrepancy. 

The good news is that the differences between 
Wilson and Clover spectroscopy on the dynamical 
ensemble are entirely consistent with the p value 
of about 200 MeV that we found in the quenched 
study. In all observables we studied, the relative 
change upon inclusion of the Clover term was less 
than or approximately 10%, i.e. (pa). In fact, in 
comparing chiral intercepts and slopes of the vec- 
tor, nucleon, and A, we were only able to resolve 
a difference between the Clover and Wilson val- 
ues of the vector intercept, i.e. the mass of the 
p. This can be seen by examining figures [|, |^, 
and |6[ in which Wilson and Clover chiral plots of 
the masses of these three states as a function of 
pseudoscalar mass are superposed. The p plot is 
the only one in which there is a clear separation 
between the Clover and Wilson data. More good 
news is that we found consistency with experi- 
ment at the 5-10% level in the light hadron spec- 
trum, since the Clover a _1 's we obtained from dif- 
ferent light observables ranged from 1.8-2.1 GeV. 



aM v vs. aM p , z 




o.o 0.2 0.4 0.6 

aM PS 2 



Figure 4. Chiral plot of vector mass on dynami- 
cal ensemble comparing Wilson(o) and Clover(o) 
discretizations 

The bad news is that these inverse lattice spac- 
ings are still inconsistent with T spectroscopy; 



aMg vs. aM pg 




0.0 0.2 0.4 0.6 0.B 

aM PS 2 



Figure 5. Chiral plot of nucleon mass on dynam- 
ical ensemble comparing Wilson(o) and Clover(o) 
discretizations 



the 20% discrepancy has not been removed by the 
Clover discretization. Since it is not a discretiza- 
tion effect, this discrepancy is probably due to 
some other systematic error. The most likely cul- 
prit is finite volume effects; the spatial extent 
of our lattice is only 16, which corresponds to 
a 1.6-1.3 fm box depending on which lattice spac- 
ing one uses. This is much larger than bottomo- 
nium, but probably too small for light hadrons. 
Squeezing the light hadronic states would raise 
their masses, resulting in spuriously low inverse 
lattice spacings. Systematic errors in the chiral 
extrapolation could also be causing the problem 
- adding a term linear in M„ (which might be 
induced by partial quenching effects) to the chi- 
ral ansatz can significantly reduce the extrapo- 
lated mass values. Another possibility is that the 
dynamical fermions are not light enough to act 
as true up and down quarks. Whatever the rea- 
son for the discrepancy, higher quality dynami- 
cal calculations, especially at larger volumes, are 
needed. It would be a tremendous accomplish- 
ment if the community could demonstrate that 
lattice QCD can simultaneously reproduce the 
spectra of both heavy and light hadron systems. 
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aM A vg. a M 
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Figure 6. Chiral plot of A mass on dynami- 
cal ensemble comparing Wilson(o) and Clover(o) 
discretizations 



6. CONCLUSIONS 

We have shown that a Clover term, with the co- 
efficient chosen using tree-level tadpole-improved 
perturbation theory, eliminates the bulk of the 
linear discretization errors in the Wilson action. 
At any value of (3, tadpole-improved Clover cal- 
culations will be much closer to the continuum 
value than the corresponding Wilson calculation. 

We have shown that the /it's for light hadron 
spectroscopy using the Wilson or Clover actions 
seem to be roughly 200-300 MeV, with the vector 
meson discretization errors slightly larger than 
those of the baryons. This means that inverse 
lattice spacings no harder than 1.3 GeV should 
be required to reduce discretization errors to the 
5% level when using the Clover action. This de- 
termination of n can be used by other groups to 
estimate the discretization errors they will face 
at a particular value of (3 and can also be used 
to determine the lattice spacings where further 
improvements of the Clover action should be at- 
tempted. 

Finally, we found that it was very difficult to 
numerically rule out linear scaling behavior when 
looking at Clover data in isolation. Only simul- 
taneous fits to both Wilson and Clover data were 



able to conclusively rule out a dominant linear be- 
havior in the Clover data in the range of lattice 
spacings studied. These fits indicated that the co- 
efficient of any sub-dominant 0(a) behavior (due 
to higher order or non-perturbative corrections to 
the perturbative Clover coefficient) is small. 
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